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Summary. This review consists of three parts: (a) what every atomic physicist 
needs to know about the physics of light nuclei; (b) what nuclear physicists can do 
for atomic physics; (c) what atomic physicists can do for nuclear physics. A brief 
qualitative overview of the nuclear force and calculational techniques for light nuclei 
will be presented, with an emphasis on debunking myths and on recent progress 
in the field. Nuclear quantities that affect precise atomic measurements will be dis¬ 
cussed, together with their current theoretical and experimental status. The final 
topic will be a discussion of those atomic measurements that would be useful to 
nuclear physics, and nuclear calculations that would improve our understanding of 
existing atomic data. 


1 Introduction 

....numerical precision is the very soul of science . 

This quote |l| from Sir D’Arcy Wentworth Thompson, considered by many 
to be the first biomathematician, could well serve as the motto of the field of 
precise atomic measurements, since precision is the raison d'etre of this disci¬ 
pline. I have always been in awe of the number of digits of accuracy achievable 
by atomic physics in the analysis of simple atomic systems^. Nuclear physics, 
which is my primary field and interest, must usually struggle to achieve three 
digits of numerical significance, a level that atomic physics would consider a 
poor initial effort, much less a decent final result. 

The reason for the differing levels of accuracy is well known: the theory 
of atoms is QED, which allows one to calculate properties of few-electron 
systems to many significant figures ||. On the other hand, no aspect of nu¬ 
clear physics is known to that precision. For example, a significant part of the 
“fundamental” nuclear force between two nucleons must be determined phe¬ 
nomenologically by utilizing experimental information from nucleon-nucleon 
scattering^], very little of which is known to better than 1%. In contrast to 
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that level of precision, energy-level spacings in few-electron atoms can be mea¬ 
sured so precisely that nuclear properties influence significant digits in those 
energies^. Thus these experiments can be interpreted as either a measure¬ 
ment of those nuclear properties, or corrections must be applied to eliminate 
the nuclear effects so that the resulting measurement tests or measures non¬ 
nuclear properties. That is the purview of this review. 

The single most difficult aspect of a calculation for any theorist is assigning 
uncertainties to the results. This is not always necessary, but in calculating 
nuclear corrections to atomic properties it is essential to make an effort. That 
is just another way to answer the question, “What confidence do we have 
in our results?” Because it is important for atomic physicists to be able to 
judge nuclear results to some degree, this discussion has been slanted towards 
answers to two questions that should be asked by every atomic physicist. The 
first is: “What confidence should I have in the values of nuclear quantities that 
are required to analyze precise atomic experiments?” The second question 
is: “What confidence should I have that the nuclear output of my atomic 
experiment will be put to good use by nuclear physicists?” 


2 Myths of Nuclear Physics 

Every field has a collection of myths, most of them being at least partially true 
at one time. Myths propagate in time and distort the reality of the present. A 
number of these are collected below, some of which the author once believed. 
The resolution of these “beliefs” also serves as a counterpoint to the very 
substantial progress made in light-nuclear physics in the past 15 years, which 
continues unabated. 

My myth collection includes: 

• The strong interactions (and consequently the nuclear force) aren’t well 
understood, and nuclear calculations are therefore unreliable. 

• Large strong-interaction coupling constants mean that perturbation the¬ 
ory doesn’t converge, implying that there are no controlled expansions in nu¬ 
clear physics. 

• The nuclear force has no fundamental basis, implying that calculations 
are not trustworthy. 

• You cannot solve the Schrodinger equation accurately because of the 
complexity of the nuclear force. 

• Nuclear physics requires a relativistic treatment, rendering a difficult 
problem nearly intractable. 

All of these myths had some (even considerable) truth in the past, but 
today they are significant distortions of our current level of knowledge. 
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3 The Nuclear Force 

Most of the recent progress in understanding the nuclear force is based on 
a symmetry of QCD, which is believed to be the underlying theory of the 
strong interactions (or an excellent approximation to it). It is generally the 
case that our understanding of any branch of physics is based on a framework 
of symmetry principles. QCD has “natural” degrees of freedom (quarks and 
gluons) in terms of which the theory has a simple representation. The (strong) 
chiral symmetry of QCD results when the quark masses vanish, and is a 
more complicated analogue of the chiral symmetry that results in QED when 
the electron mass vanishes. The latter symmetry explains, for example, why 
(massless or high-energy) electron scattering from a spherical (i.e., spinless) 
nucleus vanishes in the backward direction. 

The problem with this attractive picture is that it does not involve the 
degrees of freedom most relevant to experiments in nuclear physics: nucle¬ 
ons and pions. It is nevertheless possible to “map” QCD (expressed in terms 
of quarks and gluons) into an “equivalent” or surrogate theory expressed in 
terms of nucleon and pion degrees of freedom. This surrogate works effectively 
only at low energy. The small-quark-mass symmetry limit becomes a small- 
pion-mass symmetry limit. In general this (slightly) broken-symmetry theory 
has TOttC 2 <C A, where the pion mass is m T c 2 = 140 MeV and A ~ 1 GeV is 
the mass scale of QCD bound states (heavy mesons, nucleon resonances, etc.). 
The seminal work on this surrogate theory, now called chiral perturbation the¬ 
ory (or yPT), was performed by Steve WeinbergQ, and many applications 
to nuclear physics were pioneered by his student, Bira van Kolck[Q. From 
my perspective they demonstrated two things that made an immediate im¬ 
pact on my understanding of nuclear physicsQ: (1) There is an alternative 
to perturbation theory in coupling constants, called “power counting,” that 
converges geometrically like (Q/A) N , where Q ~ m T c 2 is a relevant nuclear 
energy scale, and the exponent N is constrained to have N > 0; (2) nuclear 
physics mechanisms are severely constrained by the chiral symmetry. These 
results provide nuclear physics with a well-founded rationale for calculation. 

This scheme divides the nuclear-force regime in a natural way into a long- 
range part (which implies a low energy, Q, for the nucleons) and a short-range 
part (corresponding to high energy, Q, between nucleons). This is indicated 
in Fig. (1), which is a cartoon of the potential between two nucleons meant 
only to indicate significant regions and mechanisms. Since yPT is effective 
only at low energies, we expect that only the long-range part of the nuclear 
force can be treated successfully by utilizing only the pion degrees of freedom. 
This would be the region with r > b. We need to resort to phenomenology 
(i.e., fitting to nucleon-nucleon scattering data) to treat systematically the 
short-range part of the interaction (r < b). 

The long-range nuclear force is calculated in much the same way that 
atomic physics calculates the interactions in an atom using QED. Both are 
illustrated in Fig. (2). The dominant interaction between two nucleons is the 
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Fig. 1. Cartoon of the nuclear potential, V(r), showing regions of importance 


exchange of a single pion illustrated in Fig. (2b) (One-Pion-Exchange Poten¬ 
tial or “OPEP”) and denoted 14 . Its atomic analogue is one-photon exchange 
in Fig. (2a) (containing the dominant Coulomb force). Because it is such an im¬ 
portant part of the nuclear potential, it is fair to call 14 the “Coulomb force” 
of nuclear physics. Smaller contributions arise from the two-pion-exchange 
potential in Fig. (2e) (called “TPEP”), which is the analogue of two-photon 
exchange between charged particles shown in Fig. (2d). There is even an ana¬ 
logue of the atomic polarization force in Fig. (2g), where two electrons simulta¬ 
neously polarize their nucleus using their electric fields. The nuclear analogue 
involving three nucleons simultaneously is illustrated in Figs. (2h) and (2i), 
and is called a three-nucleon force ||. Although relatively weak compared to 
14 (a few percent), three-nucleon forces play an important role in fine-tuning 
nuclear energy levels. The final ingredient is an important short-range inter¬ 
action (which must be determined by phenomenology) shown in Fig. (2c) that 
has no direct analogue in the physics of light atoms. Just as one can exchange 
three photons, three-pion-exchange is possible and is depicted in Fig. (2f). 

It is worth recalling that the uncertainty principle tells us that exchanging 
light particles produces longer-range forces and exchanging heavier particles 
produces shorter-range forces. Thus OPEP has a longer range than TPEP, as 
illustrated in Fig. (1). Many mechanisms have been proposed for the short- 
range part of the nuclear force, such as heavy-meson exchange, for example. 
Any meson-exchange mechanism produces singular forces, which are regular¬ 
ized to make them finite. However one chooses to do this, the part of the 
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nuclear force inside b must be adjusted to fit the nucleon-nucleon scattering 
data, and no individual parameterization of the short-range force is intrinsi¬ 
cally superior (i.e., it doesn’t matter how you do it). 



Fig. 2. First- and second-order (in a, the fine structure constant) atomic interactions 
resulting from photon exchange are shown in the left-most column, where solid 
lines are electrons, wiggly lines are photons, and the shaded line is a nucleus. The 
analogous nuclear interactions resulting from pion exchange are shown in the middle 
column, where solid lines are nucleons and dashed lines are pions. Nuclear processes 
involving short-range interactions (shaded vertical areas) are shown in the right-most 
column, together with a three-pion-exchange interaction 


How all of this works in practice is indicated in Fig. (3). Imagine that 
you throw away all of the nuclear potential inside r = b (with b chosen to be 
1.4 fm) in Fig. (1), keeping only the tail of the force between two nucleons. 
Now compute a phase shift (the 3 Po, for example). This very modest physics 
input predicts the basic shape of the phase shift (dashed line) as a function 
of energy. This variation with energy is a consequence of the small pion mass 
(compared to the energy scale in the figure). What is missing in this curve is a 
smooth (negative) short-range contribution that grows roughly in proportion 
to the energy. We can fill in the missing short-range interaction inside r = b by 
adding a potential term specified by one short-range parameter. This produces 
the dotted line, which is a rather good fit, and adding two more terms (solid 
line) produces a nearly perfect fit to the experimental results. Fixing the 
short-range part of the potential looks very much like making an effective- 
range expansion. All useful physics is specified by a few parameters, and the 
details are completely unimportant. 
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Fig. 3. 3 P 0 phase shift (in degrees) calculated only with the OPEP tail for r > b 
(dashed line), and with one (dotted) and three (solid) short-range-interaction terms 
added. The experimental results are indicated by separate points with error barsft(i[ 


What are the consequences of exchanging a pion rather than a photon? 
The pseudoscalar nature of the pion mandates its spin-dependent coupling 
to a nucleon, and this leads to a dominant tensor force between two nucle¬ 
ons. Except for its radial dependence, the form of V n mimics the interaction 
between two magnetic dipoles, as seen in the Breit interaction, for example. 
Thus we have in nuclear physics a situation that is the converse of the atomic 
case: a dominant tensor force and a smaller central force. In order to grasp 
the difficulties that nuclear physicists face, imagine that you are an atomic 
physicist in a universe where magnetic (not electric) forces are dominant, and 
where QED can be solved only for long-range forces and you must resort to 
phenomenology to generate the short-range part of the force between electrons 
and nuclei. 

Although this may sound hopeless, it is merely difficult. The key to han¬ 
dling complexities is adequate computing power, and that became routinely 
available only in the late 1980s or early 1990s. Since then there has been ex¬ 
plosive development in our understanding of light nuclei. Underlying all of 
these developments is an improved understanding of the nuclear force. It is 
convenient to divide nuclear forces and their history into three distinct time 
periods. 
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First-generation nuclear forces were developed prior to 1993. They all con¬ 
tained the one-pion-exchange force, but everything else was relatively crude. 
The fits to the nucleon-nucleon scattering data (needed to parameterize the 
short-range part of that force) were indifferent. 

Second-generation forces were developed beginning in 1993 Q. They were 
more sophisticated and generally very well fit to the scattering data. As an 
example of how well the fitting worked, the Nijmegen group (which pioneered 
this sophisticated procedure) allowed the pion mass to vary in the Yukawa 
function defining 14, and then fit that mass. They also allowed different masses 
for the neutral and charged pions that were being exchanged and found JuJ 

m n ± = 139.4(10) MeV, (1) 

m*. o = 135.6(13) MeV, (2) 

both results agreeing with free pion masses (m ff ± = 139.57018(35) MeV and 
m „o = 134.9766(6) MeV |l2|]). It is both heartening and a bit amazing that 
the masses of the pions can be determined to better than 1% using data taken 
in reactions that have no free pions! This result is the best quantitative proof 
of the importance of pion degrees of freedom in nuclear physics. 

Third-generation nuclear forces are currently under development. These 
forces are quite sophisticated and incorporate two-pion exchange, as well as 
14. All of the pion-exchange forces (including three-nucleon forces) are being 
generated in accordance with the rules of chiral perturbation theory. One 
expects even better fits to the scattering data. This is clearly work in progress, 
but preliminary calculations and versions have already appeared [p~3|. 


4 Calculations of Light Nuclei 

Having a nuclear force is not very useful unless one can calculate nuclear 
properties with it. Such calculations are quite difficult. Until the middle 1980s 
only the two-nucleon problem had been solved with numerical errors smaller 
than 1%. At that time the three-nucleon systems 3 H and 3 He were accurately 
calculated using a variety of first-generation nuclear-force models Q. Soon 
thereafter the a-particle ( 4 He) was calculated by Joe Carlson, who pioneered 
a technique that has revolutionized our understanding of light nuclei: Green’s 
Function Monte Carlo (GFMC)[pA|. 

The difficulty in solving the Schrodinger equation for nuclei is easily un¬ 
derstood, although it was not initially obvious. Nuclei are best described in 
terms of nucleon degrees of freedom. Nucleons come in two types, protons and 
neutrons, which have nearly the same masses and can be considered as the up 
and down components of an “isospin” degree of freedom. If one also includes 
its spin, a single nucleon thus has four internal degrees of freedom. Two nu¬ 
cleons consequently have 16 internal degrees of freedom, which is roughly the 
number of components in the nucleon-nucleon force (coupling spin, isospin 
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Fig. 4. GFMC calculations of the binding energies of the levels (labelled by their 
spins and parities) of light nuclei with as many as ten nucleons. These calculations 
use a common Hamiltonian and have a numerical uncertainty on the order of 1%. 
Heavy shaded lines to the left are calculated energies (with errors), while light shaded 
lines to the right are experimental energies. The label “IL2” refers to the Illinois-2 
model of the three-nucleon force that is used in all of the calculations together with 
the Argonne Vis two-nucleon force 
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and orbital motion in a very complicated way). To handle this complexity one 
again requires fast computers, and that is a fairly recent development. 

The GFMC technique has been used to solve for all of the bound (and 
some unbound) states of nuclei with up to 10 nucleons. One member of this 
collaboration (Steve Pieper|l^|) calculated that the ten-nucleon Schrodinger 
equation requires the solution of more than 200,000 coupled second-order 
partial-differential equations in 27 continuous variables, and this can be ac¬ 
complished with numerical errors on the order of 1%! A subset of the results 
of this impressive calculation are shown in Fig. (4)|l9|. 

Although the nucleon-nucleon scattering data alone can predict the bind¬ 
ing energy of the deuteron ( 2 H) to within about 1/2%, the experimental bind¬ 
ing energy is used as input data in fitting the nucleon-nucleon potential. The 
nuclei 3 H and 3 He (not shown) are slightly underbound without a three- 
nucleon force, and that force can be adjusted to remedy the underbinding. 
This highlights both the dominant nature of the nucleon-nucleon force and 
the relative smallness of three-nucleon forces, which is nevertheless appropri¬ 
ate in size to account for the small discrepancies that result from using only 
nucleon-nucleon forces in calculations of nuclei with more than two nucleons. 

Once the 3 H binding energy is fixed, the binding energy of 4 He is then accu¬ 
rately predicted to within about 1%. The five-nucleon systems (not shown) are 
unbound, but their properties are rather well reproduced. The six-nucleon sys¬ 
tems are also well predicted. There are small problems with more neutron-rich 
nuclei (compare 9 Li with ' Li or 8 He with 6 He or 4 He), but only 3 adjustable 
parameters in the three-nucleon force allow several dozen energy levels to be 
quite well reproducedBecause nuclei are weakly bound systems, there are 
large cancellations between the (large) potential and (large) kinetic energies, 
leaving small binding energies. The results shown in Fig. (4) are quite re¬ 
markable, especially given that small (fractional) discrepancies in the energy 
components lead to large effects on the binding energies. 

We note finally that power counting can be used to show that light nuclei 
are basically non-relativistic, and relativistic corrections are on the order of a 
few percent. Power counting is a powerful qualitative technique for determin¬ 
ing the relative importance of various mechanisms in nuclear physics. 


5 What Nuclear Physics Can Do for Atomic Physics 

With our recently implemented computational skills we in nuclear physics 
can calculate many properties of light nuclei with fairly good accuracy. This 
is especially true for the deuteron, which is almost unbound and is com¬ 
putationally simple. Although nuclear experiments don’t have the intrinsic 
accuracy of atomic experiments, many nuclear quantities that are relevant to 
precise atomic experiments can also be measured using nuclear techniques, 
and usually with fairly good accuracy. 
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What quantities are we talking about? The nuclear length scale is set 
by R ~ 1 fm = 10~ 5 A. The much larger atomic length scale of ao ~ 1 A 
means that an expansion in powers of R/ciq makes great sense, and a typical 
wavelength for an atomic electron is so large compared to the nuclear size that 
only moments of the nuclear observables come into play. This also corresponds 
to an expansion in a, the fine-structure constant, and m e R , where m e is the 
electron mass. This is a rapidly converging series. 

For processes that have nuclear states inside loops (such as polarizabilities) 
the excitation energies of those states play a significant role. Although states 
of any energy can be excited in principle, in practice the effective energy 
of (virtual) excitation for light nuclei (call it u>n) is within a factor of two 
of 10 MeV, except for the more tightly bound a-particle, which also has a 
smaller radius as a consequence. This number follows from the uncertainty 
principle and the fact that nucleons in a light nucleus have a radius of about 
2 fm. The deuteron’s weak binding generates the lowest values, which is about 
6 MeV for the deuteron’s electric polarizability («e ~ 1 /u>n)- Using the value 
of ui n = 10 MeV, we find u>nR ~ 1/10, which is a reasonably small expansion 
parameter. 


Table 1 . Orders in a where various contributions to the Lamb shift for S-states 
have been calculated. The label “f.s.” denotes a contribution from nuclear finite 
size or nuclear structure. Once nuclear physics enters a process at a given order, 
higher orders will also have nuclear corrections. A indicates that although a 
complete calculation of nuclear contributions has not been made, such contributions 
are expected. Names refer to the person who first calculated the leading-order term 
of that type. References and the meanings of other labels are given in the text 


Process 

2 4 

a a 

a 5 

a 6 

NR Coulomb 

Bohr f.s. 

f.s. 

f.s. 

Rel. Coulomb 

Dirac 


f.s. 

Recoil 

Darwin 

- 

- 

Nucl. Structure 


f.s. 

- 

Vacuum Pol. 


Ueliling 

f.s. 

Radiative 


Bethe 

f.s. 


At what levels do various nuclear mechanisms affect the Lamb shift? 
The (lowest) orders in a that receive nuclear contributions (for S-states) 
are sketched in Table (1). The various mechanisms are divided into static 
Coulomb (both non-relativistic and relativistic), recoil (inverse powers of the 
nuclear mass, M), nuclear structure, vacuum polarization, and radiative pro¬ 
cesses. The nuclear effects are conveniently divided into two categories: those 
that directly involve only the properties of the nuclear ground state, and those 
that involve virtual excited states and are traditionally called “nuclear struc¬ 
ture.” A radius is a good example of the former, while a polarizability is the 
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prototype of the latter. Calculational techniques are quite different for these 
two categories. 

It is beyond the scope of this review to list detailed formulae and extensive 
references to past work. I strongly recommend the recent review of || , which is 
extremely well organized. An entire section is devoted to nuclear contributions, 
and these are listed in their Table (10) with references and numerical values 
for the hydrogen atom. A sketch of how these contributions scale is given 
below together with some of the more recent references. 

The leading-order non-relativistic energy is simply the Bohr energy of 
order a 2 . Nuclear finite-size contributions of non-relativistic type (i.e., gener- 
ated by the Schrodinger equation) begin for S-states in order (Za) 4 |£G| and 
are proportional to R 2 ; they have also been calculated in order ( Za ) 5 and 
(Za) 6 || @. The Dirac energy has a leading-order (Za) 4 term, while the 
nuclear finite-size contributions of relativistic type begin in order (Za) 6 and 
are proportional to R 2 . A recent calculation exists for deuterium |23||. The 
non-relativistic finite-size corrections of order (Za) 5 and (Za) 6 are tiny for 
electronic atoms (they contain higher powers of m e R, which is very small), 
but are not necessarily small for muonic atoms (m^R is about 1 for most 
light nuclei), which was the original motivation for developing them. P-state 
finite-size effects begin in order (Za) 6 and are of both relativistic (~ R 2 ) and 
non-relativistic ( R 4 ) types. 

The most important nuclear-structure mechanism is the electric polariz¬ 
ability (which has a long history and will be discussed in more detail later), 
and this generates a leading contribution of order a 2 (Za) 3 . Coulomb correc¬ 
tions of order a 2 (Za) 4 were developed in the context of a greatly simplified 
model of the polarizability in muonic atoms j24| (which would not be applica¬ 
ble to electronic atoms). 

The Uehling mechanism for vacuum polarization is of order a(Za) 4 , while 
the first nuclear corrections are of order a(Za) 5 [|2C], 27, ^8[ ^9|. The leading- 
order radiative process is also of order a(Za) 4 , while the nuclear finite-size 
corrections begin in order a(Za) 5 [^8| |2fJ. Both of these nuclear corrections 
are proportional to R 2 . Recoil corrections have a long and interesting history 
that predates the Schrodinger equation (C. G. Darwin derived the leading 
term of order (Za) 4 using Bohr-Sommerfeld quantization; see the references 
in J3(J). To the best of my knowledge no published calculation exists for the 
nuclear-finite-size recoil corrections, which begin in order (Za) 5 , although the 
techniques of Q lead to a result proportional to R 2 /M, which should be very 
small. We note finally the hadronic vacuum polarization, which (although not 
nuclear in origin) is generated by the strong interactions [fl2| . 

One quantity through which nuclear size manifests itself is the nuclear 
charge form factor (the Fourier transform of the nuclear ground-state charge 
density, g), which is given by 
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F( q) = j d 3 r q(y) exp(iqr) 91 Z(l-^-(r 2 ) ch +-■-)-± q a qP Q a(3 +-■ ■ , (3) 

where q is the momentum transferred from an electron to the nucleus, Q a P 
is the nuclear quadrupole-moment tensor, Z is the total nuclear charge, and 
(r 2 ) c h is the mean-square radius of the nuclear charge density. These moments 
should dominate the nuclear corrections to atomic energy levels because |q| 
in an atom is set by the (very small) atomic scales. Using F to construct the 
electron-nucleus Coulomb interaction, one obtains 


Vc(r) 


Za 

r 


2irZa 

3 


<r 2 )ch£ 3 (r) 


Qa (3 (S • r) 2 — S 2 ) 
2r 3 5(25 - 1) 


(4) 


where S is the nuclear spin operator and Q is the nuclear quadrupole moment 
(which vanishes unless the nucleus has spin S > 1). The Fourier transform of 
the nuclear ground-state current density has a similar expansion 


J (q) = J d 3 r J (r) exp(iq ■ r) = -iq x [i (1 - ^-{ r2 )M + '' 0 + ' ■ • > (5) 

where fi is the nuclear magnetic moment and (r 2 )M is the mean-square ra¬ 
dius of the magnetization density. The first term generates the usual atomic 
hyperhne interaction. 


Table 2. Values of the root-mean-square charge and magnetic radii and the 
quadrupole moment (if nonvanishing) of the nucleons and various light nuclei ob¬ 
tained by nuclear experiments, together with a selected reference. If two values are 


given, the second value 

is 

Tat obtained by an atomic or 

molecular measurement 

Nucleus 

{r 2 )H 2 (fm) 

ref. 

( r2 >M 2 ( fm ) 

ref. 

Q (fm 2 ) ref. 

H 

0.880(15) 
0.883 (14) 


m 

1 

0.836 (9) 

0 

— 

2 H 

2.130(10) 

0 

2.072 (18) 

0 

0.282 (19) B 

0.2860 (15) pq, B 

3 h 

1.755 (87) 



1.84(18) 

I 

37 

1 - 

3 He 

1.959 (34) 
1.954 (8) 


41] 

1.97(15) 


37 

- 

4 He 

1.676 (8) 


45)] 

- 




- 

Nucleon 

(r 2 ) ch (fm 2 ) 

ref. 

( r2 )M 2 ( fm ) 

ref. 


n 

-0.1140 (26) 4S] 

0.873 (11) 

144] 
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Electron-nucleus scattering experiments are the primary technique used 
to measure moments of nuclear charge and current densities that are relevant 
to atomic physics |3^|, and some appropriate values of these quantities are 
tabulated in Table (2). An exception is the measurement of the deuteron’s 
quadrupole moment (Q = 0.282(19) fm 2 ) obtained by scattering polarized 
deuterons from a high-Z nuclear target at low energy Q. This result is con¬ 
sistent with the molecular determination (Q = 0.2860(15) fm 2 ) [fi9[ |~i7| , but its 
error is an order of magnitude larger. Although there is no reason to believe 
that the (tensor) electric polarizability of the deuteronQ plays a significant 
role in the H-D (molecular) quadrupole-hyperfine splitting that was used to 
determine Q, that correction was not included in the analysis. It was included 
in the analysis of the nuclear measurement. 

I highly recommend the recent review of electron-nucleus scattering by 
Ingo Sick[[17|, which contains values of the charge and magnetic radii of light 
nuclei. That review not only lists the best and most recent values of quantities 
of interest, but discusses reliability and technical details for those who are 
interested. One result from that review is listed in Table (2) and is important 
for the discussion below. The errors of the tritium ( 3 H) radii are nearly an 
order of magnitude larger than those of deuterium. Of all the light nuclei 
tritium is the most poorly known experimentally, although the charge radius 
can now be calculated with reasonable accuracy. 



Fig. 5. The direct two-photon process is shown in (a), the crossed-photon process in 
(b), and “seagull” contributions in (c). The seagulls reflect non-nucleonic processes 
and terms necessary for gauge invariance. In these graphs the double lines represent 
a nucleus, the single lines an electron, the wiggly lines a (virtual) photon. The 
shading represents the set of nuclear excited states. The loop momentum is q, and 
integrating over this momentum sets the scales of the nuclear part of the process 


In addition to moments of the nuclear charge and current densities, various 
components and moments of the nuclear Compton amplitude can play a sig¬ 
nificant role. Mechanisms that contribute to the polarizabilities are shown in 
Fig. (5). The direct (sequential) exchange of photons and the crossed-photon 
process are shown in (a) and (b), while the “seagull” process is shown in (c). 
The latter mechanism is required by gauge invariance in any model of hadrons 
with structure. The exchange of pions between nucleons generates such terms, 
for exampleQ. Because these are loop diagrams, they involve an integral 
over all momenta (q), and this sets the nuclear scales of the problem. The nu- 
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clear size scale (J?), the electron mass ( m e ), and the average virtual-excitation 
energy (tDjv, appropriate to the shaded part of the line in (a) and (b) that in¬ 
dicates excited nuclear states) determine the generalized polarizabilities p7[ . 
The process is dominated by the usual electric and magnetic polarizabilities 
and their logarithmic modifications^]. 

Specific examples are the (scalar) electric polarizability, cte, and the nu¬ 
clear spin-dependent polarizability (~ S). The latter interacts with the elec¬ 
tron spin to produce a contribution to the electron-nucleus hyperfine splitting. 
There exists a recent calculation of the latter for deuterium ][l9f . Values of the 
nuclear electric polarizability for light nuclei obtained from calculations or 
experiments are listed in Table (3). There are either calculations or measure¬ 
ments ofa E for 2 Hg| III ll), 3H H and 3 HeH|liJ,||, ^ 4He ® & 

With the exception of the 3 He experimental results, there is reasonable con¬ 
sistency. The decreasing size of the polarizabilities for heavier nuclei is caused 
by their increased binding. The a-particle has more than 10 times the binding 
energy of the deuteron, and its polarizability is an order of magnitude smaller. 


Table 3. Values of the electric polarizability of light nuclei, both theoretical and 
experimental, where the latter have been determined by nuclear exp eriments. No 
uncertainties were given for the 3 H, 3 He, and 4 He calculations in |52l l53|, but they 
are likely to be smaller than about 10%. The 4 He result was used in analyses of 
muonic He El, Ell 


Nucleus al? 10 (fm 3 ) ref. 


' (fm 3 


ref. 


H 


0.6328 (17) 


0.61 (4) 
0.70 (5) 


[5C| 
[ 511 


3 H 


0.139 


He 


0.145 


0.250 (40) 
0.130(13) 


12L 

[ 35 ] 


4 He 


0.076 


[531 0.072(4) 


[ 24 | 


The physics of hyperfine splittings is in general rather different from the 
physics that contributes to the Lamb shift. It should therefore be no surprise 
that the nuclear physics that contributes to hyperfine splittings is also quite 
different; it is also more complicated than its Lamb-shift counterpart. The 
dominant nuclear physics that we discussed previously was the physics of the 
nuclear charge density, in the form of moments of the static charge density 
(i.e., radii) and electric dipole moments that contribute to polarizabilities. 

The primary nuclear mechanism in hyperfine splittings is the magnetic 
interaction caused by the nuclear magnetization density. This density is not 
as well understood as the charge density. The primary reason is that the same 
mesons whose exchange binds nuclei together also contribute to the nuclear 
currents if they carry a charge. The pions that we discussed earlier generate a 
very important component of that current [Q. The reason for the dichotomy 
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between nuclear charges and currents can be understood by imagining that 
charged-meson exchange between nucleons is instantaneous. In this limit we 
know that the transmitted charge is always on a nucleon. In other words only 
nucleon degrees of freedom matter, which is the normal situation for the charge 
density. The power counting that we discussed earlier states that corrections 
to the nuclear charge operator are small (~ 1%), and include a type that 
vanishes for instantaneous meson exchanges. That is not the case for the cur¬ 
rent, however, since any flow of charge (even from a virtual meson) produces 
a current that is not simply related to nucleon degrees of freedom, and that 
current can couple to photons. These meson-exchange currents (often denoted 
“MEC”) can be as large as those generated by the usual nuclear convection 
current. Various tricks can be used to eliminate part of our ignorance, but 
the nuclear current density is less well understood than the nuclear charge 
density. Atomic hyperfine splittings provide us with an excellent opportunity 
to learn about nuclear currents in a very different setting. 

Although most of the hyperfine experiments in light atoms were performed 
decades ago, there has recently been renewed theoretical interest, and the ac¬ 
curacy of the QED calculations is sufficient to extract nuclear information^]. 
The differences between the QED calculations and the experimental results 
can be interpreted as nuclear corrections, and those are significant, as indi¬ 
cated in Table (4). The S-state results in this table (presented as a ratio) have 
been taken from Table (1) of the recent work of Ivanov and Karshenboim |57|| . 


Table 4. Difference between hyperfine experiments and QED hyperfine calculations 
for the nth S-state of light hydrogenic atoms times n 3 , expressed as parts per million 
of the Fermi energy. This difference is interpreted as nuclear contributions to the 
hyperfine splitting ]57|]. A negative entry indicates that the theoretical prediction 
without nuclear corrections is too large 



n\Ei - - 

rQED 

-^hfs 

)/E f (ppm) 


State 

H 

2 h 

3 H 

3 He + 

IS 

-33 

138 

-38 

222 

2S 

-33 

137 

- 

221 


Hyperfine structure is generated by short-range interactions. The dom¬ 
inant Fermi contribution (Ep) arises from a (5-function, and that produces 
a dependence on the square of the electron’s nth S-state wave function at 
the origin, |</ n (0)| 2 , which is proportional to 1/n 3 . Most nuclear effects have 
the same dependence (~ 1/n 3 ), which has been removed from the results in 
Table (4). The IS and 2S results are seen to be consistent at this level of 
accuracy, with IS experimental results typically being much more accurate. 

More calculations of the nuclear contributions to hyperfine splittings in 
light atoms are badly needed if we are to use this information to learn about 
the currents in light nuclei. These contributions come in the form of Zemacli 
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moments |38| (ground-state quantities) and spin-dependent polarizabilities 
(discussed above). There exists a considerable literature on the latter sub¬ 
ject dating back 50 years. The recent work of Mil’shtein and Khriplovich|l9|] 
has pointed out a serious defect in that older work. Although the leading-order 
terms are essentially non-relativistic in origin (for the nucleons in a nucleus), 
the sub-leading-order terms are not, and require relativity (for the nucleons) in 
order to obtain a correct result. This is not terribly surprising, since the same 
physics that enters that polarizability also enters the Gerasimov-Drell-Hearn 
sum rule (5^], which requires relativity at the nucleon level |60|), and is a topic 
of considerable current interest in nuclear physics [^TJ. The calculations of [^9| 
suggest that deuterium at least can be understood using fairly simple nuclear 
models. This needs to be checked using more sophisticated models. We note 
that the Zemach correction Js^] adds to the ratio in Table (4), improving the 
agreement between experiment and theory for H and 3 H. The large positive 
value of that ratio for deuterium suggests a large polarizability correction, 
which is confirmed by J49|. 


6 The Proton Size 


One recurring problem in the hydrogen Lamb shift is the appropriate value of 
the mean-square radius of the proton, (r 2 ) p , to use in calculations. Some older 
determinations162 disagree strongly with more recent ones |63| . As shown in 
(3), the slope of the charge form factor (with respect to q 2 ) at q 2 = 0 deter¬ 
mines that quantity. The form factor is measured by scattering electrons from 
the proton at various energies and scattering angles. 

There are (at least) four problems associated with analyzing the charge- 
form-factor data to obtain the proton size. The first is that the counting 
rates in such an experiment are proportional to the flux of electrons times the 
number of protons in the target seen by each electron. That product must be 


measured. In other words the measured form factor for low q 2 is (a — + - • •)> 

where b/a = ( r 2 ) p . The measured normalization a (not exactly equal to 1) 
clearly influences the value and error of (r 2 ) p . Most analyses unfortunately 
don’t take the normalization fully into account, and [^| estimates that a 
proper treatment of the normalization of available data could increase (r 2 )\/ 2 
by about 0.015 fm and increase its error, as well. In an atom, of course, the 
normalization is precisely computable. 

Another source of error is neglecting higher-order corrections in a (i.e., 
Coulomb corrections), and Q demonstrates that this increases (r 2 ) p ^ 2 by 
about 0.010 fm. A similar problem in analyzing deuterium data was resolved in 
[BGl. Another difficulty that existed in the past was a lack of high-quality low- 
cf data. The final problem is that one must use a sufficiently flexible fitting 
function to represent F( q), or the errors in the radius will be unrealistically 
low. All of the older analyses had one or more of these flaws. 
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Most of the recent analyses |33[ |33|, [m| are compatible if the appropriate 
corrections are made. An analysis by Rosenfelder |33] contains all of the ap¬ 
propriate ingredients, and he obtains (r 2 )j/ 2 = 0.880(15) fm. There is a PSI 
experiment now underway to measure the Lamb shift in muonic hydrogen, 
which would produce the definitive result for (r 2 ) p ^0|. One expects the 

results of that experiment to be compatible with Rosenfelder’s result. Extrac¬ 
tion of the proton radius p5| from the electronic Lamb shift is now somewhat 
uncertain because of controversy involving the two-loop diagrams. These dia¬ 
grams are significantly less important in muonic hydrogen, where the relative 
roles of the vacuum polarization and radiative processes are reversed. 


7 What Atomic Physics Can Do for Nuclear Physics 

The single most valuable gift by atomic physics to the nuclear physics commu¬ 
nity would be the accurate determination of the proton mean-square radius: 
(r 2 ) p . This quantity is important to nuclear theorists who wish to compare 
their nuclear wave function calculations with measured mean-square radii. In 
order for an external source of electric field (such as a passing electron) to 
probe a nucleus, it is first necessary to “grab” the proton’s intrinsic charge 
distribution, which then maps out the mean-square radius of the proton prob¬ 
ability distribution in the wave function: {r 2 ) wfn- Thus the measured mean- 
square radius of a nucleus, (r 2 ), has the following components: 

(r 2 ) = (r 2 ) w f n + (r 2 ) p + r 2 )„ + r 2 )..., (6) 

where the intrinsic contribution of the N neutrons has been included as well 
as that of the Z protons, and (r 2 )... is the contribution of everything else, in¬ 
cluding the very interesting (to nuclear physicists) contributions from strong- 
interaction mechanisms and relativity in the nuclear charge density |67^] . Be¬ 
cause the neutron looks very much like a positively charged core surrounded 
by a negatively charged cloud, its mean-square radius has the opposite sign 
to that of the proton, whose core is surrounded by a positively charged cloud. 
It should be clear from (6) that (r 2 ) p (which is much larger than (r 2 ) n ) is 
an important part of the overall mean-square radius. Its present uncertainty 
degrades our ability to test the wave functions of light nuclei. 

The next most important measurements are isotope shifts in light atoms or 
ions. Since isotope shifts measure differences in frequencies for fixed nuclear 
charge Z, the effect of the protons’ intrinsic size cancels in the difference. 
This is particularly important given the current lack of a precise value for the 
proton’s radius. The neutrons’ effect is relatively small and can be rather easily 
eliminated, and thus one is directly comparing differences in wave functions, 
or of small contributions from (r 2 ).... Isotope shifts are therefore especially 
“theorist-friendly” measurements, since they are closest to measuring what 
nuclear theorists actually calculate. 
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Fig. 6. Cartoon of the 2 H - H isotope shift, illustrating how the effect of the finite 
size of the proton (shaded small circle) in deuterium is cancelled in the measurement. 
The finite size of the neutron (open small circle) and the electromagnetic interaction 
mediated by the strong-interaction (binding) mechanism (indicated by the jagged 
line between the nucleons) do affect the deuteron’s charge radius (see text) 


Precise isotope-shift measurements have been performed for 4 He - 3 He[411 
and for 2 H - 1 H (D-H)||. A measurement of 6 He - 4 He is being undertaken!681 
at ANL. Gordon Drake has written about and strongly advocated such mea¬ 
surements in the Li isotopes(69|. These are all highly desirable measurements. 
Because the 3 H (tritium) charge radius currently has large errors, in my opin¬ 
ion the single most valuable measurement to be undertaken for nuclear physics 
purposes would be the tritium-hydrogen ( 3 H - 4 H) isotope shift. An extensive 
series of calculations using first-generation nuclear forces found (r' 2 )Hft for 
tritium to be 1.582(8) fm, where the “error” is a subjective estimate [|7Cl[ . This 
number could likely be improved by using second-generation nuclear forces, 
although it will never be as accurate as the corresponding deuteron value, 
which we discuss next. 

The D-H isotope shift in the 2S-1S transition reported by the Garching 
group [|| was 


Av = 670 994 334.64(15) kHz. (7) 

T T T T 

Most of this effect is due to the different masses of the two isotopes (and 
begins in the first significant figure, indicated by an arrow). The precision is 
nevertheless sufficiently high that the mean-square-radius effect in the sixth 
significant figure (second arrow) is much larger than the error. The electric 
polarizability of the deuteron influences the eighth significant figure, while the 
deuteron’s magnetic susceptibility contributes to the tenth significant figure. 
It becomes difficult to trust the interpretation of the nuclear physics at about 
the 1 kHz level, so improving this measurement probably wouldn’t lead to an 
improved understanding of the nuclear physics. 
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Analyzing this isotope shift and interpreting the residue (after applying 
all QED corrections) in terms of the deuteron’s radius leads to the results J7l[ 
in Table (5). The very small binding energy of the cleuteron produces a long 
wave function tail outside the nuclear potential (interpretable as a proton 
cloud around the nuclear center of mass), which in turn leads to an easy and 
very accurate calculation of the mean-square radius of the (square of the) wave 
function. Subtracting this theoretical radius from the experimental deuteron 
radius (corrected for the neutron’s size) determines the effect of (r 2 )... on the 
radius. Although this difference is quite small, it is nevertheless significant and 
half the size of the error in the corresponding electron-scattering measurement 
(see Table (2)). The high-precision analysis in Table (5) of the content of the 
deuteron’s charge radius would have been impossible without the precision of 
the atomic D-H isotope-shift measurement. This measurement has given nu¬ 
clear physics unique insight into small mechanisms that are at present poorly 
understood Q . 


Table 5. Theoretical and experimental deuteron radii for pointlike nucleons. The 
deuteron wave function radius corresponding to second-generation nuclear potentials 
and the experimental point-nucleon charge radius of the deuteron (i.e., with the 
neutron charge radius removed) are shown in the first two columns, followed by the 
difference of experimental and theoretical results. The difference of the experimental 
radius with and without the neutron’s size is given last for comparison purposes Jf3[ 


« (fin) 

ex P (r 2 )p{ 2 (fm) 

difference (fm) 

A(r 2 )i /2 (fm) 

1.9687(18) 

1.9753(10) 

0.0066(21) 

-0.0291(7) 


8 Summary and Conclusions 

Nuclear forces and nuclear calculations in light nuclei are under control in a 
way never before attained. This progress has been possible because of the great 
increase in computing power in recent years. Many of the nuclear quantities 
that contribute to atomic measurements have been calculated or measured to a 
reasonable level of accuracy, a level that is improving with time. Isotope shifts 
are valuable contributions to nuclear physics knowledge, and are especially 
useful to theorists who are interested in testing the quality of their wave func¬ 
tions for light nuclei. In special cases such as deuterium these measurements 
provide the only insight into the size of small contributions to the electro¬ 
magnetic interaction that are generated by the underlying strong-interaction 
mechanisms. In my opinion the tritium-hydrogen isotope shift would be the 
most useful measurement of that type. One especially hopes that the ongoing 
PSI experiment is successful in measuring the proton size via the Lamb shift 
in muonic hydrogen. 
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